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Infrared interferometric demonstrations with the University of California, Berke- 
ley’s infrared spatial interferometer (ISI) on Mt. Wilson explore the potential of in- 
frared and optical astrometry for deep space tracking, reference frame development, 
and DSN science. Astrometric data taken and analyzed over the last 5 years from the 
ISI have shown that instrumental and atmospheric effects limit current demonstra- 
tions. The benefits of sensitivity upgrades, which were performed in 1991 and 1992, 
have been demonstrated by comparing point-to-point phase fluctuations for the fall 
1989 and fall 1992 observing epochs. This comparison showed that point-to-pomt 
phase fluctuations due to tropospheric and quantum noise, for optimal integration 
times of 0.2 sec, are approaching the 0.1-cycle level needed to reliably connect the 
interferometric phase. The increase in sensitivity, coupled with that arising from 
very recent hardware upgrades, will greatly enhance phase-connection capabilities 
necessary for astrometry in the presence of atmospheric refractivity fluctuations. 
The current data set suggests that atmospheric fluctuations on Mt. Wilson during 
the best seeing are dominated by a low-lying component, approximately 25 m high, 
which may be minimized with in situ calibration in the future. During poor seeing 
conditions that currently prohibit the interferometric phase connection necessary for 
astrometry, fluctuations seem to be generated by atmospheric inhomogeneities at 
much higher altitudes above Mt. Wilson. Data taken over the last year suggest that 
the ISI will soon be able to achieve 50- to 100-nrad astrometry in a single observing 
session, employing current ground-based laser distance interferometer calibrations 
to minimize atmospheric effects. 


I. Introduction 

Infrared and optical astrometry support three areas 
of potential interest to the Deep Space Network (DSN): 
tracking of laser-carrying spacecraft; reference frame de- 


velopment, including locating solar system objects in in- 
frared/optical, radio, and planetary reference frames; and 
DSN astrometric science of the future. Because optical 
telemetry is being considered for high data-rate trans- 
mission within the next few decades [1], development of 
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an infrared or optical metric tracking capability will re- 
duce or eliminate the need for bulky radio equipment 
on laser-carrying spacecraft. Even in the absence of op- 
tical telemetry, the development of infrared and optical 
reference frames, and the location of solar system ob- 
jects in those frames, will enable more accurate target- 
relative tracking with radio systems utilizing radio-optical 
frame ties. Optical measurements of Global Position- 
ing System (GPS) satellites constitute one example of a 
radio-optical frame- tie technique. Such measurements are 
planned for the Table Mountain Ronchi telescope in the 
near future. Current DSN science frequently involves the 
measurement of electromagnetic phase or amplitude at ra- 
dio wavelengths. Similar measurements at infrared and 
optical wavelengths — such as solar deflection experiments, 
free of charged-particle effects, and asteroid imaging — will 
probably be performed by the DSN of the future. 

Stellar astrometric demonstrations with the University 
of California, Berkeley’s infrared spatial interferometer 
(ISI) bear on each of the three areas of DSN interest men- 
tioned above. ISI astrometric demonstrations, which help 
to determine the limiting accuracy of ground-based in- 
frared or optical astrometry, provide a valuable assessment 
of the limiting accuracy of tracking, reference frame work, 
and science at infrared and optical wavelengths. For exam- 
ple, atmospheric refractivity fluctuations studied with ISI 
are likely to be the dominant astrometric error for track- 
ing spacecraft, galactic, or solar system targets. Refrac- 
tivity fluctuations are similar in spatiotemporal character- 
istics at infrared and optical wavelengths. In the presence 
of refractivity fluctuations, interferometry at the longer 
infrared wavelengths enables larger apertures than those 
practical at optical wavelengths. On the other hand, two- 
color methods at optical wavelengths may help to reduce 
atmospheric effects [2]. Data from the ISI have already 
helped to characterize the Mt. Wilson atmosphere, as will 
be discussed in Section III. Descriptions of Mt. Wilson at- 
mospheric behavior may well apply to other mountain-top 
astronomical sites. 

In addition to the description of atmospheric fluc- 
tuation errors, tracking the trajectories of asteroids is 
another potential application of infrared interferometry 
in the DSN. Asteroid astrometry benefits both reference 
frame development and science [3]. Figure 1 shows the 
approximate flux of an asteroid located at 1 astronomical 
unit (AU) from the Sun. Various simplifications regarding 
the angular distribution of thermal flux from the asteroid 
make the figure approximate at the 30-percent level, but 
the illustrated wavelength dependence of the radiation is 
much better than that. The figure applies to the two cases 


noted, a 10-km-diam asteroid at 0.05 AU from the Earth or 
a 200-km-diam asteroid at 1 AU from the Earth. The ISI 
operates at 11 jivn and, with recent sensitivity upgrades, 
should be able to make astrometric measurements of ob- 
jects with flux in the 1000- Jy range. Optical or infrared 
measurements of GPS satellites will help to locate the as- 
teroids in the radio VLBI frame used for angular spacecraft 
tracking. Astrometry of asteroid trajectories will also help 
to determine their masses [3]. 

Ground-based demonstrations with the ISI will suggest 
improvements to the infrared and optical astrometric tech- 
niques, as well as the nature of future space-based tracking 
systems. This article describes recent LSI data acquisition 
and analysis that help to establish an instrumental and 
atmospheric error budget. In the next section, the current 
instrumental characteristics of the ISI and the restrictions 
they impose on astrometric performance will be discussed. 
Section III shows characteristics of the Mt. Wilson atmo- 
sphere that must be addressed in future demonstrations, 
and Section IV describes the potential for upgrading both 
the instrumentation and analysis algorithms to achieve as- 
trometry below the 50-nrad level. 


II. Instrumental Sensitivity and Systematics 

A single telescope of the ISI is shown in Fig. 2. Each of 
the two telescopes consists of a steerable flat mirror (on the 
right), a focusing parabolic mirror, optics, and electronics 
that detect the infrared signal and convert it to a 2-GIIz- 
wide radio signal for cross-correlation. Helium-neon laser 
distance interferometers (HeNe LDUs) in each telescope 
monitor the path lengths between the flat and parabolic 
mirrors and behind the flat mirror to the optics table. The 
sensitivity of the ISI for infrared astrometry is limited by 
the noise in the output of the heterodyne receivers, which 
is about twice the quantum limit, or equivalent to a single- 
sideband system temperature of approximately 2600 K. An 
ISI interferometric time series for the star Alpha Orionis is 
shown in Fig. 3. These data were taken on September 8, 
1992, on a 13-m baseline. There are points in the time 
series when a cycle may have been incorrectly assigned, 
but by and large the phase seems properly connected. 

A general requirement on the sensitivity of astromet- 
ric devices is that the white (quantum or thermal) noise 
be low enough to allow short integration times as com- 
pared to those characteristic of appreciable propagation 
media effects. For interferometric devices, both white 
noise and low-frequency-peaked propagation noise [4,5] 
contributions to the phase must be small enough to allow 


2 


phase connection. Usually rms phase variations <0.1 cy- 
cle guarantee reliable phase connection. These qualitative 
statements are illustrated quantitatively in Fig. 4, which 
shows the point-to-point rms phase variation as a function 
of integration time for data taken from Fig. 3 and from Al- 
pha Orionis on October 12, 1989. Between the two epochs 
shown in Fig. 4, a number of sensitivity enhancements were 
made, including the installation of an improved hetero- 
dyne detector in one telescope of the ISI. Low-frequency- 
peaked atmospheric effects increase as the integration in- 
terval increases. Quantum noise from the heterodyne de 
tectors, on the other hand, is white or frequency inde- 
pendent and should therefore average and produce smaller 
rms phase variations with increased integration time. The 
decreasing trend with increasing integration times until 
about 0.2 sec is presumably due to quantum noise, while 
the increasing trend from 0.2 sec on is associated with 
atmospheric refractivity fluctuations. The minimum rms 
phase fluctuation — i.e., the point at which the system noise 
and atmospheric trends cross— determines the optimal in- 
tegration time. As can be seen from the figure, sensitivitity 
upgrades improve performance by lowering the point-to- 
point fluctuation associated with the optimal integration 
time between about 0.2 and 0.13 cycle. If there were no 
rise in Fig. 4 due to atmospheric fluctuations, then in- 
tegration times could be extended indefinitely and poor 
sensitivity could be tolerated. During nights of excellent 
seeing, for example on October 10, 1989, integration times 
of 1 sec or more were used to connect phase. Figure 4 thus 
illustrates the relationship between sensitivity and atmo- 
spheric fluctuations. Based on Fig. 4, “good seeing” will 
be defined from here on as atmospheric phase fluctuations 
less than 0.1 infrared cycle on 0.2-sec time scales. Reliable 
phase connection requires slightly better sensitivity than 
that evidenced in Fig. 4 for typical Mt. Wilson observing 
conditions. Future hardware upgrades, including installa- 
tion of a new heterodyne detector in the other telescope of 
the ISI, should improve the potential for phase connection 
and therefore for infrared astrometry. 


In addition to ISI sensitivity, systematic instrumental 
effects were studied by examining HeNe LDI path delay 
time series. They revealed resonances that may originate 
from the movement of the telescopes while tracking, or 
from resonances in the power supply. The resonances were 
at about 7 Hz, but sometimes multiple resonances between 
1 and 10 Hz were found. The amplitudes of those reso- 
nances are such that they just barely contribute to infrared 
phase instability, but because they occur at frequencies 
which might cause additional problems with phase con- 
nection, we will attempt to identify and minimize them 
when analyzing data taken this past summer and fall. 


III. Atmospheric Limitations in Current ISI 
Astrometric Demonstrations 

The discussion of atmospheric errors in this section as- 
sumes that phase-connection problems will be solved by 
improved ISI sensitivity, as mentioned above. W ith reli- 
able phase connection, the atmosphere still limits astro- 
metric accuracy by causing different path-length changes 
between the two telescopes of the ISI. This section dis- 
cusses the optimal use of HeNe LDI data to minimize at- 
mospheric effects and the altitude dependence of the fluc- 
tuations as derived from model calculations. 

In Fig. 5, the Alpha Orionis data from September 
8, 1992, are shown again, plotted with the telescope- 
differenced HeNe LDI path lengths, scaled to infrared cy- 
cles, and multiplied by three to demonstrate the approxi- 
mate 0.6 correlation. This correlation suggests that simple 
subtraction of the telescope-differenced HeNe LDI path 
delays from the interferometric delays is suboptimal [4]. 
The time series in Fig. 5 prompted a model calculation 
to explore the optimal utilization of the HeNe LDI data 
in the astrometric analysis and to assess the resulting as- 
trometric accuracy. A calculation of the expected level of 
correlation between HeNe LDI and interferometric path 
delays was performed using the formalism of [5], and im- 
plemented with the numerical techniques of [6]. The total 
interferometric path length due to nonzero tropospheric 
refractivity r trop (0, 4>, t) at elevation angle 9, azimuth 4> 
relative to the orientation of the ISI trailers, and time t is 
given by 


T trop (9 ,<M) = 7 'olm(^)^.0 + 2rjf e JVe(0 (0 


where r atm (0, <j>,t) is the contribution to the interferomet- 
ric path delay from the differences in atmospheric refrac- 
tivity along the electromagnetic paths from the observed 
object to each of the two telescopes of the ISI, and T He Ne ( 0 
is the ground-based, one-way, telescope-differenced HeNe 
LDI path delay due to nonzero refractivity. The HeNe 
LDI path lengths aTe defined to lie along the x-axis at 
{9,<j>) = (0,0). The factor of 2 in Eq. (1) accounts for 
the double traversal of the HeNe LDI path by the infrared 
interferometric signal: once from the flat mirror to the 
parabola (propagating to the left in Fig. 2) and once from 
the parabola through the cat’s eye to the optics table in 
back of the flat mirror (propagating to the right in Fig. 2). 

The atmospheric and HeNe LDI path delay terms in 
Eq. (1) are 
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W0,<M) - cgin6) jf [x(^2,atm(^i 0) 

X0'l, a(m j *$! 2, t))]dz 
X r 7 

T HeNe(t ) =~ [x(r2,HeNe{x)) ~ X(n,HeNe(x))]dx 
L JO 

( 2 ) 

where x(^t,afm) i s tbe refractivity at the point denoted by 
the vector rj >a j m along the line of sight at height z above 
the ith telescope. For the HeNe LDI delay, x(^i,HeNe(x)) 
is the refractivity at the vector position a distance x along 
the HeNe path for the ith telescope. In Eq. (2), h is the 
height of the turbulent atmosphere, / is the length of the 
HeNe LDI path (which is 5 m), and c is the speed of light 
in a vacuum. Using these expressions, the correlation p 
plotted in Fig. 6 is 

p = — Tiro P T HeNe > ^ 

\J < T trop X 7 HeNe > 

where <> means ensemble average. The abscissa of Fig. 6 
is h in Eq. (2). The ensemble averages of refractivity were 
evaluated using Kolmogorov -Taylor structure functions as 
in [5], with a structure constant of 4 x 10 -7 ra -i / 3 , an 
atmospheric height of 25 m, a wind speed of 1 m/sec, and 
a saturation scale of 10 m. These parameters were chosen 
because they produce temporal structure functions similar 
to those of the data of Fig. 5. The line-of-sight coordinates 
(M) = (37.6, 71 deg) were taken from the data of Fig. 5, 
and the horizontal line in Fig. 6 is p derived from the 
data of Fig. 5, assuming that a temporal average of p over 
a single scan is equal to an ensemble average over many 
scans (ergodicity). From Fig. 6, a 25-m height is inferred 
for the turbulent atmosphere. This is a model-dependent 
result, and the sensitivity of the result to departures from 
the Kolmogorov-Taylor assumptions will be discussed in 
future articles. 

The above atmospheric modelling can be used to con- 
struct an optimal least-squares estimator for the interfer- 
ometric delay at the middle of an observation interval, 
which is the fundamental quantity of astrometric interest. 
The observed interferometric delay is modelled as 

T i — 7o($! to) -j- (t — to)r(9, <f>, i 0 ) 

+ T trop(S ,<?M) (4) 


where r 0 (^o) is the delay at the reference time t Q) and the 
changes in 9 and <f> due to sidereal tracking have been ig- 
nored. The To delay and the linear delay rate f include 
contributions due to baseline and celestial source coordi- 
nates that differ from those used in the lobe rotator model 
of the ISI correlator. An optima! estimator for tq can 
be formed [7], and its error standard deviation calculated 
using the covariance of r ivov between all times t( and t-j , 
which, suppressing the 9 and <j> arguments, is given by 

cov (nrop(ti),r tro p(t J )) = <Tt r op{ti)T irop (tj)> 

~ <T t rop(ti)>< 7~trop(tj ) > 

( 5 ) 

The two ensemble averages in the second term on the 
right of Eq. (5) are nearly zero for an interferometer with 
both telescopes at the same site. Again, the formalism of 
[5] can be used to evaluate the ensemble averages that re- 
sult after inserting expressions from Eq. (2) into the first 
term on the right of Eq. (5). Figure 7 shows the calculated, 
troposphere-limited, error standard deviation for angular 
astrometry, which is c / B times the standard deviation for 
r 0 for a 13-m baseline length B as a function of scan in- 
tegration time. The upper curve shows the effect of the 
refractivity fluctuations in the absence of HeNe LDI cali- 
bration, and the lower curve shows the improvement if the 
current HeNe LDI calibrations are optimally used. For 
1000-sec scans, an approximate 100-nrad accuracy seems 
attainable for interferometry with optimal HeNe LDI cali- 
bration. For azimuths along the HeNe LDI path, the HeNe 
LDI-calibrated accuracy is improved by about 20 percent. 

It should be noted that estimation procedures, which 
are calculationally simpler that the optimal procedure de- 
scribed above, may be used in actual data analysis; a neg- 
ligible loss in astrometric accuracy may result. The for- 
mulation using Eqs. 1—5 was presented to give insight into 
the turbulent atmospheric distance scales, the nature of 
p, and the reduction in astrometric error using HeNe LDI 
calibration. 

In addition to the optimal analysis of the HeNe LDI and 
interferometric path delays, the results of Fig. 6 suggest 
that local measurements of refractivity in the first 25 m of 
the atmosphere may yield better than 100-nrad astrometry 
for observations when the HeNe LDI correlation is high. 
These local measurements could be meteorological or could 
consist of additional HeNe LDI’s that sample the vertical 
paths above the ISI. It is very important to note that the 
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above approaches, which exploit the high correlation be- 
tween the HeNe LDI and atmospheric path delays, may 
have limited utility. On many nights in the fall of 1992, 
with poorer seeing than that of September 8, 1992, HeNe 
LDI fluctuation levels not much different from those of 
Fig. 5 were observed, while the interferometric fluctuations 
were much larger than those of Fig. 5. Because the large 
fluctuations prevented reliable interferometric phase con- 
nection on those nights, the width of the power spectrum 
of the fringe amplitudes coming from the ISI correlator was 
used as the measure of the fluctuation level. The fact that 
the interferometer signal was correlated with the HeNe 
LDI delays on nights of good seeing and much less so on 
nights of poorer seeing suggests the following picture char- 
acterizing the Mt. Wilson atmosphere: During relatively 
good seeing, the atmospheric fluctuations are fairly low to 
the ground (within the first 25 m) and optimal incorpo- 
ration of HeNe LDI data and/or other ground-based cal- 
ibration strategies may yield 100-nrad-or-better infrared 
astrometry. During poor seeing, atmospheric fluctuations 
occut much higher than 25 m above the ISI and neither 
HeNe LDI data nor ground-based calibration schemes will 
be of much help. In that case, laser guide star technology 
[8] may be of use. This hypothesis is consistent with the 
picture of the atmosphere in [4], in which larger lateral 
saturation scales are attributed to nights of poorer seeing. 
The validity of this description of the atmosphere above 
Mt. Wilson and the ultimate astrometric accuracy of the 
ISI will be explored with data taken in the summer and 
fall of 1993. 


IV. Summary and Future Directions 

Applications of infrared astrometry to DSN tracking in- 
clude research into techniques for tracking infrared or op- 
tical space-borne lasers, reference frame development, and 
astrometric science at infrared or optical wavelengths. The 
commonality of atmospheric problems at infrared and op- 
tical wavelengths and the capability of larger apertures at 
infrared wavelengths make infrared interferometry a good 
tool for studying atmospheric astrometric limitations. As- 
teroid astrometry and the measurement of gravitational 
deflection close to the Sun are examples of reference frame 
development and science that may be enabled by infrared 
interferometry in the DSN. 

By comparing astrometric data from 1989 and 1992, we 
have verified that instrumental upgrades have indeed im- 
proved interferometric phase determination at short time 
scales. The point-to-point phase scatter on 0.2-sec time 
scales is about 0.13 infrared cycle for good seeing con- 
ditions. A factor of 2 improvement in the point-to-point 


phase scatter would greatly increase the reliability of inter- 
ferometric phase connection. This factor of 2 may be real- 
ized with additional hardware sensitivity upgrades. Reso- 
nances have also been identified in the instrumental HeNe 
LDI calibration path lengths. These resonances are mainly 
at 7 Hz and are just barely strong enough to affect the in- 
frared astrometric phase. 

Analysis and modelling of ISI data taken in the fall 
of 1992 suggest a two-component model for the turbulent 
atmosphere above Mt. Wilson. For good seeing condi- 
tions, a 25-m turbulent atmospheric height has been in- 
ferred, based on the correlation between HeNe LDI and 
interferometric path-length fluctuations. A small satura- 
tion scale of 10 m was also inferred from the ISI data for 
good seeing. For poorer seeing conditions, the HeNe LDI- 
interferometric correlation is weaker, suggesting that the 
turbulent atmosphere has substantial components above 
25 m. This modelling was also used to determine that the 
atmosphere-limited astrometric accuracy on an ISI 13-m 
baseline during good seeing, with optimal HeNe LDI cali- 
bration, was about 100 nrad for a 1000-sec scan. Accuracy 
at this level would be a factor of 4 better than previous 
infrared astrometric results [9]. This level of astrometric 
accuracy has yet to be demonstrated on multiple sources 
with the ISI. Single-source phase traces and atmospheric 
modelling have been used to infer the potential accuracy 
of the ISI. 

Reliable phase connection, optimal application of HeNe 
LDI calibration, and the resulting 100-nrad single-source 
astrometry will be demonstrated with data taken in the 
next observing season. Multiple-source astrometry will 
also be attempted. This will require developing acquisition 
or analysis techniques to resolve cycle ambiguities between 
observations of multiple sources. Future astrometric im- 
provements may include local monitoring of atmospheric 
effects to improve accuracy on nights of good seeing. Such 
monitoring could involve a combination of meteorological 
sensors and new HeNe LDI paths. In the more distant fu- 
ture, the development of an infrared reference frame and 
asteroid astrometry will be pursued. A direct-detection (as 
opposed to heterodyne) system is being considered for the 
ISI on 5-year time scales. Direct detection would greatly 
increase the usable ISI bandwidth and sensitivity, enabling 
astrometry in a much wider variety of seeing conditions. 

It should be noted that because of the apparent small 
saturation scale during periods of good seeing, increases 
in baseline length will yield almost proportionate improve- 
ments in accuracy. Data taken on azimuths closer to that 
of the HeNe LDI will be more effectively calibrated by op- 
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timal application of HeNe LDI data. Considering these 
factors, astrometric accuracy of the order of 50 100 nrad 
is probably characteristic of future ISI single-observation 
performance in the absence of the above-mentioned me- 
teorological calibrations. Because of the small saturation 


scales, troposphere-induced errors for observations from 
the same session, or from different sessions, should be un- 
correlated. By averaging results over many observations, 
astrometric accuracies of much better than 50 100 nrad 
should be possible. 
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Fig. 1. Black-body calculation of approximate flux density for an 
asteroid located 1 AU from the Sun, as a function of wavelength. 
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